Quantum spin mixing in a binary mixture of spin-1 atomic condensates 
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We study quantum spin mixing in a binary mixture of spin-1 condensates including coherent 
interspecies mixing process, using the familiar spinor condensates of 87 Rb and 23 Na atoms in the 
ground lower hyperfine F = 1 manifolds as prototype examples. Within the single spatial mode 
approximation for each of the two spinor condensates, the mixing dynamics reduce to that of three 
coupled nonlinear pendulums with clear physical interpretations. Using suitably prepared initial 
states, it is possible to determine the interspecies singlet-pairing as well as spin-exchange interactions 
from the subsequent mixing dynamics. 
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I. INTRODUCTION 

A topical area in physics today concerns the control 
and manipulation of the spinor degrees of freedom asso- 
ciated with electrons or atoms. Two highly visible sub- 
fields attracting tremendous theoretical and experimen- 
tal interests are spintronics in condensed matter systems 
[l[ and atomic spinor quantum gases 0. The latter sys- 
tem become available due to the technical breakthrough 
of optical trapping, which provides equal confinement for 
all atomic Zeeman components of a fixed F. As a re- 
sult, spin-related phenomena are exhibited and detected 
in cold atoms, including various quantum phases [H-Il] 
and quantum magnetism studies [9(, the observation of 
spin domain formation [ToL IT]] ] , as well as the dynamics 
of spin mixing [f| , and spin squeezing fl2l [l3j , etc. 

According to the formulation of atomic spinor conden- 
sates 0-[l], the order parameter for a condensate in the 
hyperfine F state is generally described by a spinor of 
2F + 1 components, strongly influenced by their under- 
line atom-atom interactions. Within the low energy limit 
of interests to atomic quantum gases, when modeled by 
contact interactions, atom-atom interactions are required 
to be invariant with respect to both spatial and spin ro- 
tations, reflecting the nature of s-wave interactions. De- 
pending on the values of spin-dependent interaction pa- 
rameters, the ground state of a spinor condensate can be 
ferromagnetic or anti-ferromagnetic (polar) for F = 1 0- 
_while an additional cyclic phase appears when F — 2 
Higher spin cases are generally more complicated 
with limited experimental access. 

Law et al. pioneered the study of atomic spin mix- 
ing dynamics [5[ . They first adopted numerical approach 
studying quantum spin mixing in the absence of an ex- 
ternal magnetic (B-) field @. Subsequent theoretical and 
experimental efforts have contributed to the observation 
and control of the coherent quantum dynamics, otherwise 
rarely visible in many body systems [14l422l | . 

In the semiclassical picture, using mean-field approx- 
imation and adopting the single spatial mode approxi- 
mation (SMA) HHII, coherent spin mixing dynamics in 



a spin-1 condensate is described by a nonrigid pendu- 
lum, displaying periodic oscillations and resonance be- 
havior in an external B-fAeld [13, HH . This picture proves 
to be widely popular with experimentalists and led to 
many successes pdMl9j . Analogous efforts were applied 
to spin-2 condensates, for instance, in the hi gher hy per- 
fine manifold of the ground state 87 Rb atoms [191 4221 . An 
interesting application suggested by Saito et al. [26[ pro- 
vides a practical method for determining the unknown 
spin coupling parameters (polar or cyclic) relying on the 
mixing dynamics with suitably prepared initial states. 

Recently, several groups investigate intensively mix- 
tures of spinor condensates (27rl33| , whose properties are 
reasonably well understood, both when an external B- 
field is absent or present. As before for a single species 
spinor condensate, semiclassical mean field approxima- 
tions are adopted and the full quantum approach is lim- 
ited to atom number dynamics in the restricted spatial 
modes of the condensates. The ground state properties 
for the mixture, is to a large degree, determined by the 
yet unknown interspecies spin exchange interaction. If 
it is antiferromagnetic and is sufficiently strong, interest- 
ing pha ses, such as highly fragmented ground states arise 
[29 . |3(J | . Additionally, there exists the so-called broken- 
axisymmetry phase in the presence of an external B-ficld 
[3ll | . Within the degenerate internal state approxima- 
tion [34[ , which considers atomic interaction potentials as 
coming from contributions of potential curves associated 
with the coupled electronic spins of the two valence elec- 
trons: one for each type of atoms (taken as alkali atoms 
for simplicity) [13, [HI, H^|. The interspecies singlet- 
pairing interaction vanishes as all interspecies interaction 
parameters are determined by a total of only two scatter- 
ing lengths for the electronic singlet and triplet channels 
respectively. This approximation provides a zeroth order 
estimates for the 87 Rb and 23 Na atom mixture we study. 
Experiences with spin exchange interactions within each 
species show otherwise, i.e., the need for more atomic 
interaction parameters. 

We therefore propose to develop analogous spin mix- 
ing dynamics as in F — 2 spinor condensates. We 
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will calibrate the interspecies singlet-pairing interactions 
with suitably prepared initial states as in spin-2 conden- 
sates [26[. Additionally, we find that inter-species spin- 
exchange interaction can also be determined analogously. 




II. THE MODEL OF A BINARY SPIN-1 
CONDENSATE MIXTURE 

The binary mixtures of spin-1 condensates have been 
discussed in several earlier studies [2814311 ] . In addition 
to the individual Hamiltonian for each species of the two 
spinor condensates, additional contact interactions exist 
between the two species which can be decomposed into 
spin-independent and spin-dependent terms as well, de- 
scribed by Vuiri-rz) = |(a + /3Fi • F 2 +^T'o)5(r 1 - r 2 ) 
[29l l30l | with appropriate interactions parameters a, (3 
and 7 [29|, [30| . Take spin-1 condensates of 87 Rb and 
23 Na atoms as examples, the total Hamiltonian is then 
given by 
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where Hi and H 2 describe a single species system of 87 Rb 
and 23 Na atoms respectively with the interspecies inter- 
action described by H i2 . V°, Mi, pi, and gi (V 2 °, M 2 , p 2 , 
and q 2 ) respectively denote the optical trap, atomic mass, 
linear, and quadratic Zeeman shifts of a 87 Rb ( 23 Na) 
atom. Both the nuclear spins and the valence electron 
spins are the same for the two species. In the subspace 
of hypcrfinc spin angular momentum F = 1, the linear 
Zeeman shifts for both 87 Rb and 23 Na atoms are thus 
almost equal: p\ ~ p 2 (= p). ^i(r) (4>j(r)) annihilate a 
87 Rb ( 23 Na) atom at the position r. 

The F = 1 states for both 87 Rb and 23 Na atoms are 
well studied, and their respective atomic collision param- 
eters are known precisely, quote their sources of respec- 
tive ao and a 2 , which then gives ai/ 2 and fli/ 2 [37l.l38j. 
While a number of experimental and theoretical studies 
have previously addressed collisions between 87 Rb and 
23 Na atoms [HI, H6|, the most recent one by A. Pashov 
et al. [36j provides a well converged data set for sin- 
glet and triplet scattering lengths of a s = 70 (ao) and 
a f = 109(ao). This can be used to predict the required 
set of atomic intraspecies collision parameters a, (3, and 
7. What is certain concerns the value of spin exchange 
interaction 7, it will be actually strong, instead of being 
weak or vanishing. Perhaps we should consider using the 
real atomic values for some of the calculations. 




FIG. 1: (Color online). Schematic illustration of the three 
coupled nonlinear pendulums. 



We now adopt the mean-field approximation and define 
for each condensate species a mode function if}{f)/ <p{r), 
justified by the fact the spin independent density inter- 
action terms are usually much stronger than spin depen- 
dent interactions. We therefore take ^(r) = (^(f)) = 

yQp e i0 Jil){f f ) and $i(r) = ($,(r)) = \fnf ] ' e^ ^r). 
The spin dynamics arc then governed by the spin- 
dependent energy functional 
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where n lj2 = J2 3 n j ' '> m i,2 = n[ ' - nif ', m = 
0-1-26*0, r\2 = tpi + (p-i-2ip , and r/ 3 = 6-i-9i+ipi - 
tp-\. The interaction parameters are now redefined to ab- 
sorb the relevant multipliers: [3[ = pi J \ifi(f)\ 4 dr, f3' 2 = 
p 2 J\(b(r)\ 4 dr, and 08', 7') = ((3, 7 ) J \i;(r)\ 2 \cj>(r)\ 2 df. 
We note that / \ip(r)\ 2 dr = J \cf>(r)\ 2 dr = 1. When the 
two species are immersible, the overlaps between ip(r) 
and <p{r) is significantly reduced, leading to diminished 
j3' and 7', essentially reducing the system to two stand- 
alone spin-1 condensates. 

Although complicated in forms, the above Hamiltonian 
gives rise to dynamics that can be interpreted simply in 
terms of three coupled nonlinear pendulums, with three 
pairs of canonical variables: (n,Q ,T)j) and (7713 = mi — 
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W2j?73)- Their corresponding equations of motion are 
given by 
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as illustrated schematically in Fig. [TJ 



4 8£ 
V3 = t- 



(3) 



III. DETERMINING THE INTERSPECIES SPIN 
SINGLET-PAIRING INTERACTION 

When discussing spin mixing in a spin-2 condensate, 
Saito et dl. (26j proposed to determine the value of 
intra-species spin singlet-pairing interaction by choosing 
an elementary process + <-> 2 + (—2) which occurs 
only when the spin singlet-pairing interaction is non- 
vanishing. With a suitable initial state of zero magne- 
tization, the mixing dynamics is governed by coupled 
first-order ordinary differential equations, which contain 
unknown parameters like singlet-pairing interactions and 
quadratic Zeeman shifts. The analytic solutions can be 
compared with the experimental measured dynamics to 
decide the unknowns. 

Analogous approach can be taken to determine the 
value of interspecies spin singlet-pairing interaction for 
a binary mixture of spin-1 87 Rb and 23 Na atom conden- 
sates, making use of a different elementary collision pro- 
cess ^i + ^-i <-> ^-i+c&i which appears only in the 
presence of the 7' term. With an initial state 



■00 ari d re found to remain exactly zero within the 
mean-field approximation, unless dynamical instabilities 
exist. If instabilities do occur, they can be suppressed by 
tuning the quadratic Zeeman shifts qj to a large negative 
value, for instance with off-resonant microwave field (39l . 
|40| , or to a large positive value with increased uniform B- 
field. The processes <-» ^o+^o and O 

+ $0 will then be suppressed, the populations of the 
Mp = states remain at zero. In this case, the spin 
mixing dynamics of Eq. ([3]) reduce to that of a single 
pair, which takes the form, 
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and is described by a simpler energy functional 
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after neglecting a constant term — pm + (ft + ft + ft — 
7'/6)to 2 /8 + 7'ni7i2/12 + qifii + (72^2 ■ Substituting Eq. 
© into Eq. ©, we find 
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which can be integrated following the procedure of Rcf. 
[25j by solving for the equation to 3 = 0, keeping the in- 
terspecies spin-dependent interaction parameters ft and 
7' as unknown. 

From the Eq. ([5]) and assuming an initial state with 
sin 773 > 0, we infer 7' > if to 3 increases during the 
initial short time period of the spin mixing dynamics, 
and 7' < if it decreases. To determine 7', we prepare 
an initial state with 773 = n/ 2 and mi = TO2 = 0, which 
leads to £ = in the Eq. ([6]), and 
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= 



2{n\ + n 2 T V( n l ~ n l) 2 + 4C 2 n 2 ?i 2 )/(l - C 2 ). For 
C > 1, however, only two solutions —x\ and x\ exist. The 
solution for the mixing dynamics is expressed in terms of 
the Jacobian elliptic functions sn(.) and cn(.) as 



/ X27WI -C 2 x x 
TO 3 (rj = x\ sn ^£ , — J ,tor L < 1 



TO3 (t) = ii en K 



\2Ti ' x 2 

xi \ yt^/(x 2 +x 2 )(C 2 -l) 



\/xj+x 



12% 



,forC > 1, 



(9) 



where K(.) is the complete elliptic in- 
tegral of the first kind, and x 3 = 

^2(n 2 + n 2 + y/(ni - n 2 ) 2 + 4C 2 njn 2 ) / \C 2 ~^l). 

The stability of the above dynamics are confirmed with 
numerical solutions, taking the initial state as ^j/ip = 
V^(1,0,1) T /V2, $j/4>= v^(l, 0, -i) T /V2, assuming 
77i = Ti2 = n and 77 = 2 x 10 4 . We further choose ft/K. = 
-22.4893 x lCT 4 Hz and /3' 2 /h = 303.816 x 10" 4 Hz. A 
noise level at 10 -5 in the population of Mp = spin 
states of both species is also included. The B-field is 
set as large enough to suppress the intraspecies spin- 
exchange process with the quadratic Zeeman shifts sat- 
isfying qi = 40|ft|77 and q 2 = qiAEi/AE 2 , where AEi 
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FIG. 2: (Color online). Time dependent populations of each 
spin components. In the left panels of (a)-(c), the interspecies 
interaction parameters used are /?' = 5\/3[ \ and 7' = 2|/3i|. For 
the panels of (d)-(f), /3' = 5|/3(| and 7' = -2|/3(| are used, (a) 
For the 87 Rb condensate, where the solid blue line, dashed red 
line, and dotted-dash black line represent the Mf = 1,0,-1 
components, respectively, (b) As in (a), but for the 23 Na 
condensate, (c) Time dependent 7713 with solid blue line and 
red square symbols denote numerical and analytical solutions 
respectively, (d) As (a), but with 7' = — (e) As in (b), 

but with 7' = -2|/3[|. (f) As in (c), but with 7' = -2|#|. 



and AE 2 are the hyperhne splittings of 87 Rb and 23 Na 
atoms respectively. Figure [2] illustrates our numerical 
results. In Fig. I^a-c), /3' = 5\f3[\ and 7' = 2\0' 1 \ are 
used, while (3' = and 7' = — 2\j3[\ are used instead 
for Fig. [2] (d-f). The time evolution for each conden- 
sate species are shown in Fig. [5{a,d) and (b,e) respec- 
tively for 87 Rb and 23 Na atoms. The evolutions for 7773 
are shown in Fig. HJc) and (f), indeed they confirm our 
predictions based on the insights gained from analytical 
solutions that 7773 increases/decreases at the beginning 
when 7' > O/7' < 0. The numerical simulations denoted 
by solid blue lines agree well with analytical solutions 
of Eq. (JTJ]) denoted by red square symbols. We further 

note that nj' = —h_\ = —n^ = = 7713/4 with the 
initial state used in this case. As a result, we can deter- 
mine the sign of 7' from the populations of arbitrary spin 
components and species. 

Using Eq. ((9|), wc can then proceed to determine the 



value of 7' if C 
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determined from the oscillation amplitude of 777,3. Af- 
terwards, f3' becomes partially determine to within the 
following two choices 



IV. DETERMINING THE INTERSPECIES 
SPIN-EXCHANGE INTERACTION 



In the previous section, a scheme is proposed capable 
of determining the interspecies singlet-pairing interaction 
parameter 7' following spin mixing dynamics from a suit- 
ably chosen initial state. The interspecies spin-exchange 
interaction parameter /3', which is partially determined 
at the same time, will become fully determined with the 
dynamics discussed in this section. 

Equation ([2]) gives four relevant elementary spin- 
exchange processes: * + $0 *-i + *o + $0 <-> 

+ *_i+$o *o + *-i, and *i + $ «-» *o + $i, 
which can be used to determine the interspecies spin- 
exchange interaction. The hrst two processes involve 
both P' and 7' terms of the Hamiltonian in Eq. ([1}, 
while the last two processes are solely induced by spin- 
exchange interactions. 

The above four individual processes become indepen- 
dent if all other possible collision channels are suppressed. 
For example, to observe the mixing due to XPo + ^0 
"f-i + $1, ^1 = $_i = needs to be ensured at all 
times. A plausible scenario can again employ increased 
quadratic Zeeman shifts q\ and q 2 . As long as the en- 
ergy difference between the final state and the initial 
state increases, the intraspecies spin-exchange process 
*i + <-> * + *o and $1 + $_i <-» § Q + <t> are 
suppressed. They help to maintain close to zero popu- 
lations during time evolution in the corresponding spin 
state, if an initial state with <Fi = $_i=0is used. 

In the following we will describe the isolation of the 
process \&i + <&o ^0 + $1 as an example to deter- 
mine the interspecies spin-exchange interaction. An ini- 
tial state 
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is assumed, together with a sufficiently strong uniform 
external magnetic held to guarantee ^ -1 = $_i = 0. 
Since interspecies spin mixing is only induced by the 
same /3' term, and energy conservation, the spin mixing 
dynamics is then governed by the evolution of 7773 
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which can be derived from the Eq. @ , and the associated 
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energy functional 
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after neglecting a constant term — pm — + /?2 — 
/3')m 2 /8 + (/3[m + /3 2 ™2 + 9i + q 2 )m/2. Furthermore 
we can rewrite Eq. (fT2"|) as 
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The procedure to fully determine /?' goes as follows. 
First we infer the sign of j3' from the initial stage of 
the time evolution for m 3 , as in the earlier section on 
determining the sign of 7'. For an initial state with 
= V»T(1, !. 0) T /\/2 and $/<^ = 0^(1, -i, 0) T /V2, 
where (771 - 772 - J7 3 )/2 = 6»i - 6 - tpi + (fo = -t/2, 
we confirm /3' > (/?' < 0) if m 3 initially increases (de- 
creases). The actual value of /3' is determined by com- 
paring the analytic or numerical solutions using the two 
choices of (3' from the Eq. (fT0|) to experimental measure- 
ments. Again we assume n\ = n 2 = n = 2 x 10 4 , f3[/Ti = 
-22.4893 x l(T 4 Hz, (3' 2 /h = 303.816 x 10" 4 Hz, /?' = 
5\(3[\, i = 2\(3[\, Ql = 30|#|n, and q 2 = q 1 AE 2 /AE 1 , 
with the analytic solution for m 3 
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xi(x 2 — x 4 ) + {x\ — x 2 )x 4 x 2 
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sn ( c?4 — t\/d 3 /d 2 , d\) 
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where Xj=i 2,3,4 are the four roots of m 3 = 0, ar- 
ranged in desce nding order x\ > > x 2 > x 3 > x± , 
and di = y j {x\ - x 2 )(x 3 - x 4 )j (xi - x 3 )/{x 2 - x 4 ), 
d 2 = 2/y / (x 1 -x 3 )(x 2 -x i ), d 3 = [4^' 2 - 
{P[ + % + ff_ + 7 V3) 2 ]/16?i 2 , and d 4 
F(arcsin y/—(x 2 — x 4 )xi/(xi — x 2 )/x 4 , d\), with F(.) 
the elliptic integral of the first kind. 

Figure [3] show population evolutions for all spin com- 
ponents. Due to the large but unequal quadratic Zee- 
man shifts qi and q 2 , the suppression of intraspecies spin 
mixing dynamics leads to a suppressed amplitude for 
interspecies spin-exchange dynamics. As a result, the 
quadratic Zecman shifts cannot be tuned too large, oth- 
erwise they cause nonzero populations in the Mp = — 1 
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FIG. 3: (Color online). Population dynamics for all spin 
components, with the interspecies interaction parameters 
P' — 5\p[\ and 7' = 2|/3i|, and the quadratic Zeeman shifts 
qi = 30\p[ \n and q 2 = qiAE 2 /AEi. (a) For the 87 Rb conden- 
sate, where the solid blue line and dashed red line represent 
the Mf = 1,0 components, respectively, (b) As in (a), but for 
the Mf = — 1 component in dotted-dash black line. (c)/(d) 
corresponds to that in (a)/(b) respectively, but for the 23 Na 
condensate, (e) Time evolution of 7713, where solid blue line 
and red square symbols denote numerical and analytical so- 
lutions respectively. 



spin component especially for the 23 Na atoms as illus- 
trated in Fig. [3(d) . 

The other three elementary channels can also be em- 
ployed to determine the interspecies spin-exchange inter- 
action. Among them, two are capable of determining the 
combined parameter /3'— 7'/3, which can be further aided 
by a determination of the sign of j3' — 

Before conclusion, we hope to stress that the special 
mixture illustrated in this study involves a spin-1 con- 
densate with ferromagnetic interaction ( 87 Rb) and a po- 
lar spin-1 condensate ( 23 Na) with antiferromagnetic in- 
teraction. More generally the procedures we suggest for 
determine the interspecies interaction parameters remain 
applicable for mixtures with two spin-1 ferromagnetic 
condensates or two antiferromagnetic condensates. 



V. CONCLUSION 

We discuss coherent spin mixing dynamics for a bi- 
nary mixture of spin-1 condensates. Under the mean 
field approximation, the dynamics reduce to three cou- 
pled nonlinear pendulums, one for each spin-1 condensate 
as understood previously for stand-alone spin-1 conden- 
sate 



25| , and a third one for the difference in magnetiza- 
tion between the two species. By tuning quadratic Zee- 
man shifts to large enough values, they can suppress in- 
traspecies spin mixing dynamics, which results in a pure 
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interspecies spin mixing dynamics. Using suitably pre- 
pared initial states with zero populations in the Mp = 
states for both species, we can determine the value of the 
interspecies singlet-pairing interaction by comparing the 
analytic formula to experimental measurements, and at 
the same time we can partially determine the value of 
the interspecies spin-exchange interaction parameter /?'. 
Next, using an alternative initial state with zero pop- 
ulations in the Mp = —1 states of both species, and 
using the two possible values for /?' partially determined 
above, we can numerically or analytically solve the dy- 



namics and compare them with experimental results to 
determine the correct value of j3' . 
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